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with eigenvalue equation 
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of the resolvent of the Hamiltonian matrix for the system of length N with boundary 
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and Tr(T— z) -1 , together with an identity of Thouless type connecting Tr (log \T\) with the 
Hamiltonian eigenvalues for z = e 1 ^ . The results are then extended to T^T by showing that 
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§1 Introduction. 

Several discretized models are described by a Hamiltonian matrix Ti, with tridiagonal 
structure made of blocks H n = along the main diagonal, and blocks L n , L^, with 
det L n 7^ 0, respectively in the adjacent upper and lower diagonals, the blocks having size 
M x M. The diagonal matrices may describe the inner dynamics of a sequence of finite 
subsystems, and the off-diagonal matrices are the couplings among neighbouring ones. 
An important and extensively studied example is Anderson's model for electronic transport 
in a D-dimensional lattice with random potential, which for D=3 exhibits a metal-insulator 
transition [1]. The single matrices H n are random Hamiltonians for the isolated slices of 
dimension D-l and, in the simplest case, the couplings L n are proportional to the unit 
matrix, as implied by the discretization of the Laplacian. They may also be complex, in 
the presence of a magnetic field [2], or random, due to random hopping amplitudes [3]. 
Another frequently studied model is the ensemble of band random matrices [4], where H n 
is a member of GOE or GUE and L n are random and lower triangular matrices. They 
found applications in quantum chaos [5], one particle mesoscopic transport [6] and the 
propagation of two particles in disordered media [7]. The block structure also arises in the 
Fourier representation of the Floquet Hamiltonian Ht = Hq + V^e lujt + Ve~ lujt — idt, giving 
H n = Hq + nu and L = V. Block Hamiltonians have also been investigated in the context 
of matrix models [8] , the blocks being rotationally invariant for the methods to apply. 

The eigenvalue equation for TL, in block components, is: 

L n lp n+1 + H n lp n + L^i^n-l = E i>n (1-1) 



One is often interested in asymptotic properties of eigenvectors. A basic tool for this 
analysis, which exploits the recursive content of the eigenvalue equation, is the transfer 
matrix, connecting pairs of components of the vector 




The transfer matrix has size 2M x 2M and it is convenient to factorize it as follows: 



T(E) = Tjv(-E)Ejv_iTjv-i (-^)^iv-2 • • • XiT^E) (1.3) 

™>-(V V). li 

By allowing for complex values of the parameter E, for the single factors and then for the 
whole product, one obtains the important symplectic property 

T(E*y* 2 T(E) = * 2 , a 2 =(^ ~ J ) (1.4) 

where I is the unit matrix of size M. Another consequence of the factorized expression is 
detT(E) = Ilfcdet^L- 1 ], that implies | det T(E) \ = 1. 
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General theorems assert that the eigenvalues t a of transfer matrices built with random 
factors grow or decrease exponentially with N [9] , allowing the definition of characteristic 
exponents 

la {E)= Jim llog|* a (£)| (1.5) 



which, in the case of (T^T) 1 / 2 , constitute the Lyapunov spectrum of the model. The 
analytical derivation of a Lyapunov spectrum is usually extremely difficult, the alternative 
being a careful numerical work to cope with exponential instabilities. For 2x2 matrices 
(M = 1), a relevant formula by Herbert and Jones, rediscussed by Thouless and bearing 
his name, connects the single Lyapunov exponent to the eigenvalue density of the ensemble 
of Jacobi Hamiltonians [10]. However, the density is by no means a simpler problem; an 
exception is Lloyd's model, characterized by diagonal disorder with Cauchy distribution: 
in this case the analytical expression of 'j(E) is known [10]. The statistical properties of 
the Lyapunov exponent and various generalizations have been investigated extensively by 
Pendry [11]. 

One of the rare solvable examples in more than one dimension is in the work [12] by 
Isopi and Newman, who studied products of matrices all of whose entries are identically 
distributed random variables, and found analytically a "triangle law" for the Lyapunov 
spectrum; Cook and Derrida considered the case of randomly sparse matrices [13]. A 
beautiful statistical theory of transfer matrices, based on few physical contraints, has been 
introduced by Mello and others to describe transport properties in disordered multichannel 
conductors , obtaining the observed value for universal conductance fluctuations [14]. In 
general, however, the transfer matrix is a derivated object, which inherits a structure from 
the specific recurrence or dynamical equation under examination. Its statistical properties 
depend in a complicated way on fluctuating parameters that enter more naturally, for 
example, in the characterization of an ensemble of Hamiltonians. 

For the class of Hamiltonians we are considering, analytical results are lacking. Many 
extensive numerical calculations have been carried out for the Anderson model. It was by 
means of transfer matrices that Kramer and MacKinnon first exhibited the metal-insulator 
transition in 3 dimensions [15], through the different scaling behaviour in the transverse 
area M of the smallest Lyapunov exponent. A recent numerical study by Markos [16] 
provides the whole Lyapunov spectrum, which is sensitive to the transition. The Lyapunov 
spectrum of band random matrices was obtained numerically by Kottos et al. [17], with a 
discussion of finite size corrections. 

The aim of this paper is to investigate some general mathematical properties of the 
transfer matrix T(E) that arise from the block structure of a single but generic Hamiltonian 
matrix H, of size NM, corresponding to (1.1) with boundary condition i/jq = i/jn+i = 0. 

In [18] I showed that the eigenvalues of the transfer matrix are most directly related to 
those of a matrix H(z) of size NM, which in general is not Hermitian and has block 
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structure 

L{ H, 



H(z) 
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. . . Hn-i Ln-i 



(1.6) 



resulting from the eigenvalue equation (1.1) with boundary conditions specified through a 
complex parameter z: 

L N ip N+1 = zip! = -ipN (1-7) 

z 

The relation is based on the simple statement, whose proof is straightforward: 
A vector (ipi, ip2, ■ ■ ■ > V'iv) J 's an eigenvector of with eigenvalue i? if and only if 

(zipi,ipN) JS an eigenvector ofT(E) with eigenvalue z, the components ip2 • • -ipN-i be- 
ing linked to ipN and LniPn+i = zipi by (1.1). 

It implies that the characteristic polynomials of T{E) and Ti{z) are proportional, and 
eventually leads to the following " duality relation" : 

det[T(E) -z} = (-z) M det(LAr_i . . . Li)" 1 det[E - H(z)} (1.8) 

A new proof will be given in §2, after having derived an explicit representation of the 
matrix T(E) in terms of the corner blocks Gy, i,j = 1, N, of the resolvent G = (H — E)~ x 
of the Hamiltonian matrix Tt, for the system of length N. The Hamiltonian has block 
structure (1.6) with null matrices replacing the ^-dependent corner blocks, corresponding 
to the boundary requirement ipo = ipN+i = 0. 

In this paper it is noted that the derivative in the variable z gives a relation between 
resolvents 

^{nSn)-T + |logdet[ B -« W ] (1.9) 

The equation also follows from a representation of [T(E) — z\~ x in terms of the corner 
blocks of the resolvent G = (Ti.(z) — E) -1 , to be obtained in §4. Did 7i belong to an 
ensemble, it would provide access to the spectral density of T(E) by relating the difficult 
problem of averaging the resolvent of the transfer matrix, which depends on the ensemble 
parameters in a complex way, to the average of a ratio of determinants of the Hamiltonian 
itself. 

The general discussion for T(E) will be extended in §5 to cover the relevant matrix 
T(E) J[ T(E), by showing that it is itself the transfer matrix of a tridiagonal block Hamil- 
tonian, of size 2NM. 

The main results of the paper can be summarized in the following equations: relation 
(1.9) or its equivalent (3.6) among traces of resolvents, (2.3) and (4.3) that provide the 
representation of the transfer matrix in terms of corner blocks of resolvents, (3.5) and 
(5.4) that generalize Thouless' relation to M > 1 and relate Lyapunov exponents to the 
eigenvalues of the Hamiltonian matrix, the duality relation (5.6) for ther matrix T^T. 

The described results are exact and general. Hopefully, this work should provide an an- 
alytical framework for the hard task of investigating average spectral properties of transfer 
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matrices, constructed from an ensemble of Hamiltonians of this type. One more motivation 
is the recent interest in Non-Hermitian matrices, which are now extensively investigated 
mainly in the one-dimensional case M = 1 [19], precisely with the tridiagonal structure 
and boundary conditions that are here generalized. In one dimension, the eigenvalues 
are distributed along curves in the complex plane [20] and real eigenvalues correspond to 
delocalized eigenvectors [21]. The here established relationship with the transfer matrix 
makes non Hermitian matrices an obvious object for investigating localization properties 
of eigenvectors. For M = 1 this was done by Brouwer et al. [22] . 

§2 The transfer matrix and the resolvent of H. 

In this section I obtain a block representation of T(E) in terms of the corner blocks 
of the resolvent of 7i. I then prove the duality relation (1.8). 

For complex E, let us introduce the resolvent G(E) = (7i — E)~ x . It is a matrix made 
of N 2 square blocks G tJ of size M x M. The relation [G{E)i$ = G{E*) ji holds. By 
definition: 

L\_ x Gi_\j + (Hi - E)Gij + LiG i+1 j = 18^ (2.1) 

By solving the recurrence relations for j = 1 and j = N, one obtains two identities 
involving the transfer matrix 

^=T{E)( Gl A, (- I )=T(E)( G ]: N ) (2.2) 



Gn,i J V — i J ' \Gn,n / V o 

They can be joined into a matrix relation, which gives a representation of the transfer 
matrix in terms of the corner blocks of the resolvent: 

G ~+^%G ) {2M) 

One checks that the symplectic property (1.4) is identically satisfied. Note that each block 
component of T(E) is a matrix polynomial in E, and is here explicitly represented in terms 
of the resolvent of the Hamiltonian. By varying the number iV of factors in T(E), one 
obtains a collection of matrix polynomials which follow orthogonality relations that extend 
to M > 1 the familiar notion of othogonal polynomials of Jacobi matrices [23]. 
To the end of deriving the duality relation, I first prove 

detG-jv = det[L 1 ...L N _ 1 ]- 1 det[E -H] (2.4) 

Proof: a vector (ipi, . . . , i^n) is an eigenvector of 7i with eigenvalue E if and only if it solves 
(1.1) with boundary conditions V'jv+i =^0 = 0. These conditions, by eqs. (1.2) and (2.3b), 
are equivalent to the requirement = G± N ip\. By construction T(E) is a polynomial in 

E of degree iV with matrix coefficients; then detG^ is a polynomial in E of degree 
NM and leading term (-l) M E MN detfL^ . . .Lj -1 ].' The polynomials det[E - H] and 
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detfG^jy] share all zeros and are therefore proportional by a numerical factor obtained 
from comparison of the leading terms.* 

I now give a proof of the duality relation (1.8), alternative to that provided in [18]. 
Proof: by writing 7i(z) = H + V(z), where V(z) is zero everywhere except for the corner 
blocks Vi t N = \l and Vn,i = zl, one calculates: 

det[E - H{z)) = det[E - H] det [J + G(E)V(z)] = 

I + zG hN iG M \ (2-5) 



= det[E-H}det( I + zG ^ N ) 



On the other hand, by using algebraic properties of determinants, one obtains from eq. 
(2.3b): 

det[T(£) -z)=z M det[G^ N ] det ( 1 J ^ ) (2.6) 

By taking into account property (2.4), the duality relation (1.8) follows. • 
§3 The duality relations. 

I here discuss some consequences of the duality relation, 

det[T(£) -z} = (-z) M det(L 7V -i . . . £i) _1 det[E - H{z)} (3.1) 

An identity for the inverse transfer matrix is obtained from the adjoint of (3.1), with the 
aid of the properties T(£*)+ = -a 2 T(E)- 1 a 2 and H(z*Y = H(l/z): 

det[T( J E)- 1 -z} = (-z) M det(L) v _ 1 . . . lJ)" 1 det[E - H{l/z)\. (3.2) 

The product of the two identities immediately yields an identity which was used in [18] to 
study the structure of bands and the dynamics of eigenvalues of the Hermitian Hamilto- 
nians H{e iLp ): 

det[T(E)+T(E)- 1 -(z + -)} = | det(Ljv-i • • • £i)|" 2 det[E-H(z)\ det[E -H(l/z)} (3.3) 

z 

A simple general consequence of (3.1) is that, for ImE ^ 0, the transfer matrix T{E) has no 
eigenvalues on the unit circle, since the right term of (3.1) never vanishes for a Hermitian 
matrix. More generally, this is true for E not in the union of the bands Bf~, k = 1 . . . NM, 
each one being defined as the interval of the real axis spanned by the eigenvalue Ek(<p) of 
the Hermitian matrix 7i(e*^), as <p varies in [0, 27r). 

From the symplectic property (1.4) it follows that if t is an eigenvalue of T(E) with \t\ ^ 1, 
then 1/t* is an eigenvalue of T(E*). In particular, for real E, the 2M eigenvalues of T(E) 
occur in pairs t, 1/t*, unless \t\ = 1. The number v of pairs of eigenvalues on the unit 
circle coincides with the number of bands with intersection in E [18]. 
Let us denote the eigenvalues of T(E), E real, as 

t a = e ±Xa+l9a , a = 1, ...,M-v, t b = e i6 \ b=l,...,2v 
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and write the modulus of (3.1), with real E and z = e lLp , in terms of them. After some 
simple algebra: 



M—v 2u ,^ x 

Y[ (2coshA a -2cos(6' a -v?)) JJ 2 sin ( -\0 b - (p\ j = 

= \det(L N - 1 ...L 1 )\- 1 \det[E -H(e lv )]\ (3.4) 

By taking the logarithm of it and integrating the phase ip in [0, 2tt) one obtains a remarkably 
simple and interesting relation: 

M-u N-l 2tt 

V A a (£) = - Vlog|detL,| + — / d^log|det[£?-W( e i *')]| (3.5) 

o=l j = l ^ ^ 

This formula is exact, and is valid for a single matrix. In the large N limit and in a statis- 
tical context, one would have the average behaviour (A a ) = N^ a , where 7 a is independent 
of N. The right hand side of (3.5) would be evaluated by means of the ensemble and 
(^-averaged density of eigenvalues of H(e' l(p ). The equation would then provide a general- 
ization to M > 1 of the Thouless relation between the Lyapunov spectrum and the average 
eigenvalue density of the Hamiltonian ensemble [10]. 

If instead we take in (3.1) the derivative in the variable z, by using the property 
d/dz det(A + zl) = det(A + zI)tr(A + zl) -1 , we obtain equation (1.9). Also the derivative 
in the right hand side can be computed, and gives the following final formula, where Gij 
are the blocks of size M x M that partition the resolvent G(z, E) = \H{z) — E]~ l : 



f 1 \ M m ~ . 1 
\T(E) 



TV ( ^ - )=-—- T±Gi,n + ^G NA (3.6) 



The same formula follows from a stronger result, to be given in the next section. 
§4 The transfer matrix and the resolvent of TC(z). 

A representation of the resolvent of the transfer matrix can be given in terms of the 
resolvent of the non- Hermitian matrix, G(z,E) = [H(z) - E]' 1 . Note that [G(z, E)^ = 
G(l/z*, E*)ji. With the same procedure as in §2, one obtains two identities 

(^• 1 )=T(E)( 1 A' ), ( zd ^- I )=T(E)( 1 d ^) (4.1) 
\G N ,i J \-G NA -lJ \ Gjv.iv / \-(^n,nJ 

which join into the matrix relation 

zG 1N -I *G^\ /G IlW G hl \ 



Simple steps lead to the final representation: 

1 / -Gi,jv 



T{E)-z \-\G N , N \G Ntl -\l 



(4.3) 
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which, by taking the trace, provides eq. (3.6). 

Note that the corner blocks of G can be expressed in closed form in terms of the corner 
blocks of G, by means of the Lippman Schwinger equations 

Gij = Gij H — G^iGnj + zG^nGij (4.4) 



§5 The matrix T+T. 

The general results obtained so far for transfer matrices can also be applied to the 
matrix Q(E) = T(E)^T(E), which will be shown to be itself the transfer matrix of a 
Hamiltonian built out of TC. 

The matrix Q(E) has the feature of being Hermitian and positive, therefore with real and 
positive eigenvalues. It is easy to show the validity of the property: 



Q(E*)a 2 Q(E) = a 2 (5.1) 

It follows that if t is an eigenvalue of Q{E), then 1/t is eigenvalue of Q(E*). 
While considering the factorization 

T(E)^T(E) = Ti(£)tE}r 2 (£)t . . . J^ N _ 1 T N {E)^T N {E)T, N . X . . . SiT^E) 



one notes the property that T k (E)\ constructed with H k , coincides with —T k (—E*) con- 
structed with —H k . This allows to interpret Q(E) as the transfer matrix for the solution 
of the equation JC(E)^f = 0, with matrix JC(E) = 

/ Hi — E Li \ 
L\ Ho-E 
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L 



2 



Ln-i 

L n _i H N — E —I 



-I E* — H N — Lj v _ 1 

-Ln-2 



Ln-i E* — H N -i —L n _ 2 



-Li E* — H\J 



The corresponding non-Hermitian matrix K,(E, z) entering the duality relation, is obtained 
by placing the diagonal matrices z~ x I and zl in the upper right and lower left corners 
respectively. Then, the following equation holds: 

det[T(£)tT(£) - z] = {-l) NM {-z) M \det{L N -i...L 1 )\- 2 detlC{E,z) (5.2) 

A few remarks on the spectral properties of JC(E, z), which can be easily translated for the 
matrix JC(E), are useful: 
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a) 1C(E,z)i = K(E*,l/z*); 

b) ( ° D r\K{E,z)(l " P ) = -/C(£M/c) 



n -p o y v ' ; \P o 

where P is the block matrix with nonzero blocks Pi^-i+i = I of size M, i = 1 . . . N; 

c) for real z, because of a) and b), the eigenvalues of K.{E, z) come in pairs x, —x*; 

d) the value x = does not belong to the spectrum of 1C(E, z) if z is not in the real 
positive axis, because the left side in (5.2) cannot vanish. 

Let us concentrate on the case where z = e l(p and E is real; the matrix K,(E, e l(p ) is 
Hermitian, therefore its eigenvalues are real. By point c) the eigenvalues for cp = 0, it 
are symmetric; they also mark the extrema of the bands [18]: it follows that K.{E,e l<p ) 
has as many positive as many negative eigenvalues and (— \) NM det K{E, e lv ) > 0. The 
eigenvalues of Q(E) are M — fi positive pairs (ta,^" 1 ), with t a = e Xa > 1, being 2\i 
eigenvalues equal to unity. Equation (5.2) reads, in terms of the eigenvalues of Q(E): 

(2 sin |) ^ Y[ (2 cosh A a - 2 cos v?) = JJ | det L k \~ 2 (-1) NM det 1C(E, e itf ) (5.3) 

o=l k=l 

By taking the logarithm and integrating in <p, we end with a formula of Thouless type: 

N-l 1 „2tt 

VA a (P) = -2 Vlog|detL fc | + — / cVlog[(-l)^ M det/C(£,e^)] (5.4) 
k=i Zn Jo 

For real z and E real or complex it is convenient, at least for a simplification of the 
notation, to bring the matrix JC(E, z) to another form; there is much freedom since only 
the determinant of the matrix matters. Let us choose to left and right-multiply matrix 
/C(-E, z) by unitary matrices to give: 

H-ReE + U -iV-lmE ' 
iV-ImE -H + ReE + U 

where P is the same matrix of size NM defined after (5.2), U and V are block diagonal 
matrices, each of the iV diagonal blocks having size M. The only nonzero blocks are: 
Ui : i = — z~ x )I, V\^\ = + z~ l )I and Vn,n = —I- With this transformation, we 
obtain the equivalent form of the duality relation: 

det[T(E)^T(E) - z] = (-l) NM (-z) M \ det(L Ar _ 1 . . . L^" 2 det K' (E, z) (5.6) 

The matrix JC'(E, z) is Hermitian for any E in the complex plane and real z, and it has 
the advantage of containing the matrix H in the diagonal blocks, albeit with opposite sign; 
a similar structure appears in a paper by Efetov [24]. It has the following properties: 

a) !C'(E,zy = K!{E,z*) 
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b) (J. J) (J i ^ = -k'{e\i/z) 

c) (J Q _ 7 ) =/?(£*, -1/z) 
§6 Conclusions. 

In the present paper and in [18] I have considered a class of Hamiltonians character- 
ized by a block Jacobi structure which is shared in many interesting models of quantum 
disordered transport. For a single Hamiltonian matrix I have obtained exact relations 
that allow to describe spectral properties of transfer matrices through properties of the 
Hamiltonian itself. The identities involve general boundary conditions that imply a close 
connection between transfer matrices and Non-Hermitian Hamiltonians, and make the 
statistical analysis less involved, since the statistical ensemble is usually defined for the 
Hamiltonian. 

Some equations, like (1.9), are suited for the supersymmetric technique. Two problems 
arise, that have already been considered in the literature [25]: 1) the need of a special for- 
malism for the determination of the density of complex eigenvalues from the knowledge of 
the average resolvent, 2) the " Hermitianization" procedure for representing ratio of deter- 
minants, which provide by differentiation traces of resolvents, as Gaussian superintegrals. 
These problems are absent while considering the relation for T^T. 
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